Magnon Heat Transport in (Sr,La)i4Cu2404i 
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We have measured the thermal heat conductivity k of the compounds Sri4Cu24 04i and 
Ca9La5Cu2404i containing doped and undoped spin ladders, respectively. We find a huge anisotropy 
of both, the size and the temperature dependence of k which we interpret in terms of a very large 
heat conductivity due to the magnetic excitations of the one-dimensional spin ladders. This magnon 
heat conductivity decreases with increasing hole doping of the ladders. The magnon heat transport 
is analyzed theoretically using a simple kinetic model. From this analysis we determine the spin gap 
and the temperature dependent mean free path of the magnons which ranges by several thousand 
A at low temperature. The relevance of several scattering channels for the magnon transport is 
discussed. 
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Low-dimensional quantum magnets have attracted 
much attention in recent 

years among both experimental and theoretical physi- 
cists. On the one hand these compounds serve as model 
systems for a comparison between experiment and the- 
ory, since exact solutions or numerical treatments of the 
model Hamiltonians yield clear-cut predictions. On the 
other hand unusual ground states and magnetic excita- 
tions are present, in particular in quasi-one-dimensional 
spin systems. A prominent example is the quantum dis- 
ordered spin liquid state with a spin gap which is well 
established in frustrated and dimerized spin chains and 
in spin ladders 

Usually the magnetic excitations of these spin systems 
are experimentally studied by measuring spectroscopic 
or thermodynamic quantities which reveal information 
on the magnetic ground state and the excitation spec- 
tra as a function of energy and momentum. In principle, 
dispersive magnetic excitations should also contribute to 
a transport property, i.e. the thermal heat conductivity 
AC. The experimental investigation of the magnon heat 
transport could give interesting complementary informa- 
tion on the magnetic excitations, such as dissipation and 
scattering of magnons, similar as the study of electronic 
transport properties does in metals. 

Recently, several studies of thermal heat conductiv- 
ity in low dimensional spin systems have been per- 
formed and magnetic contributions to k are discussed 
in one-dimensional spin systems such as e.g. CuGeOa 
and SraCuOa as well as in two dimensional 
cuprates However, in many compounds the inter- 

pretation of the data is controversial, since an unambigu- 
ous discrimination of different contributions to k is diffi- 
cult or even impossible ]l5[ | . Very convincing experimen- 
tal evidence for a magnon heat transport has been pre- 
sented by Sologubenko et al. for (Sr,Ca)i4Cu2404i [ p^ . 
The crystal structure of this compound contains two 



quasi-one-dimensional magnetic subsystems along the c- 
axis. One subsystem is a sheet like arrangement of 
CU2O3 two-leg ladders, where the Cu 5 = ^ spins are 
strongly coupled via a Cu-O-Cu superexchange, while 
the other subsystem is an array of Cu02 S — ^ spin 
chains with weak magnetic interactions. While even the 
stoichiometric compound Sri4Cu2404i is hole doped the 
holes are predominantly located in the chains |]l7| , |l8| . 
Changing the composition i.e., the Sr, Ca, La content 
in (Sr,Ca,La)i4Cu2404i, alters concentration and distri- 
bution of the holes. This results in drastic changes of the 
magnetic properties of the chains, whereas a large spin 
gap in the ladders of the order of 400 K is observed in 
all compounds [|T9|-p2[. However, the charge transport is 
determined by the holes in the ladders which are already 
present in (Sr,Ca)i4Cu2404i. Undoped ladders are only 
found in systems containing a large amount of trivalent 
ions, e.g. La. 

In this paper we report on the thermal conductivity of 
Sri4Cu2404i and Ca9La5Cu2404i parallel and perpen- 
dicular to the chain/ladder direction. In the undoped 
ladders of Ca9La5Cu2404i the magnon contribution is 
very large and exceeds the phonon contribution by nearly 
two orders of magnitude. A simple approach to describe 
the energy transport due to magnetic excitations in the 
ladders is presented. 

We have grown single crystals of Cai4La5Cu2404i 
and Sri4Cu2404i by the traveling solvent floating zone 
method |^,|^. Using a standard steady state method 
measurements of k have been performed on pieces cut 
along the principal axes with a typical length of 2 mm 
along the measuring direction and of about 0.5 mm for 
the two other directions. The thermal gradient has been 
determined with an Au/Fe-Chromel thermocouple. 

Fig. 1 presents k of Sri4Cu2404i as a function of tem- 
perature T measured along the three crystallographic 
axes (kq, Kb, Kc). A striking anisotropy of both, ab- 
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solute value and temperature dependence of the thermal 
conductivity is apparent. Only for the 6-axis we find the 
qualitative behavior expected for phonon heat transport. 
At low T the occupation of phonon states implies an 
increase of the heat conductivity, whereas Kh decreases 
at high T due to increased phonon scattering. Ka devi- 
ates slightly from this usual phonon thermal conductiv- 
ity of insulators (see inset of Fig. 1): At high temper- 
ature we find a slight increase with increasing T. This 
small increase of Ka might be related to the complexity 
of the phonon spectrum and/or unusual scattering pro- 
cesses possibly related to the anomalous T dependence 
of the lattice constants We note that small devi- 

ations from the T dependence expected for the phonon 
heat transport as we observe for Kq in Sri4Cu2404i are 
rather common for complex transition metal oxides. It is 
still possible and meaningful to discuss this data in the 
framework of a phonon heat conductivity Kph- 
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FIG. 1. Anisotropic thermal conductivity Ka, «6 and Kc of 
Sri4Cu2404i as a function of temperature. Inset: Ka and fct 
normalized to the value at 300 K. 



This is impossible in the case of Kc, i.e. 



for 



the heat transport along the chain/ladder direction in 
Sri4Cu2404i. In this case the phonon heat transport can 
only explain the data below « 40 K where Kc exhibits a 
low T maximum very similar to and k;,. The abso- 
lute value is about 1.6 times larger than for Kq. Since 
the velocity of sound and the elastic constants do not 
show a pronounced anisotropy ||2^ one has to attribute 
the anisotropy of Kpt to anisotropic phonon scattering. 
Such an anisotropy of the phonon mean free path is also 
observed in other transition metal oxides (l^j2^ with low 
dimensional structure elements. 

The most striking observation is the behavior of Kc 
for temperatures above w 40 K. With increasing T Hc in- 
creases strongly and a second very pronounced maximum 
of the heat conductivity occurs at T ~ 140 K. Further in- 



creasing the temperature causes a sharp decrease of Kc 
and for temperatures above 250 K the heat conductiv- 
ity seems to saturate at a still rather large value of about 
18 W/Km. Both the very large values of Kc at intermedi- 
ate T and the strange temperature dependence with the 
pronounced second maximum differ drastically from the 
usual Kph- We mention that our data for Sri4Cu2404i 
agree qualitatively with those reported by Sologubenko 
et al. |]l^. There are, however, strong deviations to the 
data of Kudo et al. pq ] who report much smaller k for 
all lattice directions and T probably due to enhanced de- 
fect scattering in their samples. Comparing our data for 
Kc quantitatively with those in Ref. p6| reveals a per- 
fect agreement for the low T maximum while the second 
peak is about 30% larger in our sample. At present we 
are systematically studying different crystals in order to 
understand the origin of these quantitative deviations. 

As discussed in Ref. it is reasonable to attribute 
the high temperature maximum of Kc to magnetic excita- 
tions propagating along the spin ladders. The electronic 
heat transport as estimated from resistivity data and the 
Wiedemann-Franz-law is negligible. The T-dependence, 
anisotropy of k and the absolute value of Kc can not be 
explained by a phonon heat transport and the data for 
Cai4La5Cu2404i allow also to exclude exotic electronic 
contributions as discussed below. 

A determination and analysis of this very unusual 
magnon contribution requires a separation of a phonon 
background. Unfortunately, it is impossible to determine 
unambiguously the phonon background for Sri4Cu2404i. 
For example, at high T the absolute values and T- 
dependence are unpredictable, which is clearly demon- 
strated by the qualitatively different behavior of Ka and 
Kf). Similar problems do exist at low T, since k is deter- 
mined by strongly anisotropic phonon scattering which 
can not be determined from independent experimental 
data or from theory. This means that any determination 
of the additional magnetic contribution in Sri4Cu2404i 
leads to a large error for low T and high T, i.e. for tem- 
peratures where the total k is not much larger than Kph- 
In order to estimate the phonon background we fit the 
data at low temperatures T < 40 K with a Debye-model 
and extrapolate the behavior up to 300 K. The result of 
this description of Kph is indicated by a solid line in Fig. 1 
which is much smaller than the measured k, in the entire 
temperature range above 50 K. 

However, considering only the data for Sri4Cu2404i 
one can question any non-phononic heat conductivity at 
r>250 K. A much more reliable separation of phonon 
and magnon heat transport is possible in the case of 
Ca9La5Cu2404i whose thermal conductivities along the 
a- and c-axes are shown in Fig. 2. At low T the 
thermal heat conductivity of this compound which con- 
tains undoped spin ladders and only slightly doped spin 
chains |17 1^ is much smaller than in stoichiometric 
Sri4Cu2404i for both directions. It is straightforward 
to attribute this isotropic suppression of Kph at low T 
to structural defects. On the one hand a lattice site is 
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occupied by randomly distributed Ca^+ and La'^"'" ions 
with strongly different sizes which is known to suppress 
Kph- On the other hand x-ray diffraction studies show 
structurally disordered chains in CagLa5Cu2404i. Sur- 
prisingly, Kc strongly increases above 40 K in spite of the 
structural disorder. At intermediate and high T Kc is 
even larger than in stoichiometric Sri4Cu2404i and the 
room temperature value is comparable to that found in 
metals. In contrast to that a very small heat conductivity 
is found along the a-axis in the entire temperature range. 
Reduced size and T dependence are typical for Kph of a 
compound with many structural defects. It is apparent 
from Fig. 2 that this strong damping of Kph which is in- 
ferred from well as from the low T behavior of Kr 
enables a reliable determination of the additional contri- 
butions to Kc above « 40 K. 
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FIG. 2. Thermal conductivities of Ca9La5Cu24 04i as a 
function of temperature measured along the a- and c-axes, 
Ka and Kc, respectively. The solid line represents an estimate 
of the phonon contribution to Kc- 



mention several conclusions which can be drawn from 
the qualitative behavior of k in Ca9La5Cu2404i and, in 
particular, from the comparison to the findings in the 
stoichiometric compound. The suppression of Kph which 
coincides with an enhancement of Kc at high T gives very 
strong evidence that there are two independent 

contributions to the heat transport. Moreover, some 
thinkable origins of the large additional contributions 
to Kc in Sri4Cu2404i can be ruled out from the ad- 
ditional data. Excitations of the chains are irrelevant, 
since the magnetic and electronic properties of the heav- 
ily doped chains in Sri4Cu2404i showing charge order, 
spin gap and dimerization are completely different from 
the long range ordered chains in Ca9La5Cu2404i p9| , ^ . 
Heat transport due to collective electronic excitations 
as for example suggested for the thermal conductivity 
of the charge density wave compound K0.3M0O3 



could be relevant in the strongly doped charge order- 
ing Sri4Cu2404i whereas it is certainly irrelevant in 
Ca9La5Cu2404i due to the strongly reduced hole con- 
tent. Summarizing these arguments we have to state that 
the findings in Ca9La5Cu2404i strongly support the in- 
terpretation of the additional contribution to Kc in terms 
of a magnon heat transport along the spin ladders. This 
magnon heat conductivity Kmag is apparently not sup- 
pressed by structural disorder and, moreover, a smaller 
hole doping seems to enhance K^ag at high T. 
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FIG. 3. Magnon thermal conductivities of Ca9La5Cu2404i 
and Sri4Cu2404i as a function of temperature. 

In Fig. 3 we show Kmag for Sri4Cu2404i and 
Ca9La5Cu2404i which are derived by subtracting the 
Debye-fits of the phonon contribution from the measured 
Kc- For r< 100 K absolute value and temperature depen- 
dencies of K„iag are similar in the two compounds. How- 
ever, pronounced differences occur at higher T. The mag- 
netic heat conductivity on Sri4Cu2404i Kmag strongly 
decreases above « 150 K. This decrease is much less 
pronounced and is found at higher temperature (above 
^ 200 K) in Ca9La5Cu2404i where we find a very large 



Before we turn to a quantitative analysis of Kmag we '^•n 



100 W/Km even at room temperature. 



In order to analyze the magnon heat transport theo- 
retically we start from the simple kinetic expression 



(1) 



where Ufc, e/j, and li^ denote velocity, energy and mean 
free path of the magnetic excitations. Assuming a mo- 
mentum independent mean free path, i.e. Ik 
using the distribution function 



^mag and 



rik 



(2) 



for the triplet excitations yields the following expression 
for the magnon heat conductivity: 
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Here, A'' is number of ladders per unit area and emax is 
the band maximum of the spin excitations which is at 
approximately 200 meV in (SrT;a,La)i4Cu2404i. Note, 
that the distribution function (g) is different from a Bose 
distribution. This difference is used to account, on av- 
erage, for the hard-core constraint of no on-site double- 
occupancy for the triplet excitations. While such a form 
of the distributions function would describe the occupa- 
tion of a local triplet excitation exactly, it is only meant 
as an approximate phenomenology regarding the momen- 
tum space distribution - which is unknown for the present 
case. The primary motivation for of this type of a triplet 
distribution is to suppress unphysically large triplet den- 
sities at higher temperatures. From Eq. it is appar- 
ent that the particular form of the magnon dispersion 
does not enter the heat conductivity of a one-dimensional 
triplet gas within a kinetic description - an effect which 
has not been noticed in other studies In particular, 
the momentum or energy dependence of the magnon ve- 
locity V which emerges from Eq. |^ plays no role in our 
approach. Note that the assumption of a fc-independent 
mean free path implies a /c-dependent scattering rate 
since v depends on the wave number. While this assump- 
tion seems reasonable for defect scattering more involved 
scenarios will apply in particular to inelastic scattering. 
Within our treatment the magnon heat 

conductivity for experimentally relevant temperatures 
T <C Smax/ks mainly depends on two parameters: the 
spin gap Aiadder and the mean free path Imag- We men- 
tion, that (^) differs from an expression used by So- 
logubenko et al. for the heat-conductivity of one- 
dimensional bosons not only by the distribution func- 
tion (|^) but also by an overall factor of three accounting 
for the triplet degeneracy. 

For temperatures below 300 K Kmag does not depend 
significantly on Cmax and in the following analysis of the 
data in the framework of Eq. Jsl we use emax — 200 meV 
taken from the literature |2iy32|] . It is impossible, to de- 
termine reasonable values of the two remaining parame- 
ters Imag and l^iadder by fitting the data with the expres- 
sion in Eq. ^ without further assumptions. Therefore we 
ignore possible temperature dependencies of the spin gap 
and, moreover, assume a constant, temperature indepen- 
dent mean free path at low T, i.e. for a small number of 
thermally activated phonons and magnons. Knowing the 
spin gap from this first step of the analysis we can di- 
rectly extract the temperature dependent Imag from the 
data. 

Our assumption of a constant Imag at low T implies 
that the heat conductivity is determined by the temper- 
ature dependent activation of magnons in this tempera- 
ture range. This is in agreement with the experimental 
data which roughly follow a simple activated behavior 
as displayed in Fig. 4. Note that Eq. ^ predicts devia- 
tions from the simple activation law and taking into ac- 
count these corrections slightly improves the description 
of the data as demonstrated in Fig. 4. The spin gaps of 
^ladder = 396 K and /^ladder = 430 K wc obtain are in 



fair agreement with the results from neutron scattering. 
We conclude that both, the temperature dependence of 
Kmag and the absolute values of the spin gaps extracted 
from the data confirm our analysis. Comparing the left 
and right columns of Fig. 4 reveals that the agreement be- 
tween data and theory is much better for CagLa5Cu2404i 
than for Sri4Cu2404i. In the latter case a fair agreement 
with the theory at constant Imag can only be obtained 
in a limited temperature range 55K<T<85 K. The de- 
viations at higher temperatures, i.e. for TiJSS K signal 
a significant temperature dependence of the mean free 
path, while the differences at low temperatures T< 55 K 
are most probably due to the uncertainty of the phonon 
heat conductivity which is much larger than Kmag in this 
temperature range (see Fig. 1). The separation of Kmag 
and Kph is much more reliable in CagLa5Cu2404i and, 
indeed, the thermal conductivity for T^lOO K is well 
described by our model with constant Imag- 
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FIG. 4. Temperature dependence of the magnon ther- 
mal conductivities of Ca9La5Cu24 04i (left) and Sri4Cu2404i 
(right) at low temperatures in comparison with fit results 
to the data (Dotted lines: simple activation, solid lines: 
Id- model for n). Insets: Arrhenius representation of k, and fit 
results. 

Based on the preceding analysis of the low temperature 
behavior of Kmag it is possible to determine the tempera- 
ture dependence of the mean free path Imag (see Fig. 5). 
At low temperatures we find a very large Imag of several 
thousand A. Within the experimental uncertainty due 
to the phonon background we can not resolve differences 
of the low T saturation values of Imag for the two com- 
pounds containing spin ladders with different hole con- 
tent. There are, however, unambiguous and drastic dif- 
ferences at higher temperatures. In Sri4Cu2404i which 
contains spin ladders with a rather large hole doping the 
magnon mean free path is much smaller and we find a 
rather strong temperature dependence between 85 K and 
220 K. According to our analysis Imag saturates for higher 
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temperatures and at T = 300 K the mean free path is 
about 60 A. While the saturation of Imag above 220 K 
is clearly extractable from our data, the absolute values 
of Kmag and Irnag strongly depend on the choice of the 
background, i.e. the extrapolation of Uph, and therefore 
a much smaller Imag at 300 K is also compatible with the 
raw data in Sri4Cu2404i. For the second crystal with 
nearly undoped spin ladders a very large Imag at room 
temperature is out of question. Our analysis yields a 
magnetic mean free path larger than 500 A at 300 K. 
Moreover, the overall temperature dependent decrease 
of Imag IS less pronounccd and deviations from the con- 
stant low T value start at higher temperatures than in 

Sri4Cu2404i. 

Our results for Imag in Sri4Cu2404i are in qualita- 
tive agreement with the findings reported by Sologubenko 
et al. []l6| . There are, however, quantitative discrepan- 
cies due to the different theoretical description: for a 
given Kmag OUT analysis leads to a smaller value of Imag, 
mainly due to the more realistic distribution function for 
the magnetic excitations which we use in our approach 
(Eq. Yet, the mean free paths which we extract from 
the analysis of our data are very large. In Ref. it 
was speculated that holes are the main scatterers in the 
doped spin ladders of (Sr,Ca)i4Cu2404i. This specula- 
tion is only partially supported by our results on a crys- 
tal with undoped spin ladders. At low T our data do not 
show a significant difference between the two crystals in- 
dicating that the very large Imag at low temperature does 
not depend strongly on the hole content of the ladders. 
It is likely that the maximum magnon mean free path 
is determined by defects. Note however, that one has 
to discriminate between crystal defects and "magnetic 
defects" , since there is no correlation between the max- 
imum Imag and the maximum phonon mean free path. 
The latter is much smaller in Ca9La5Cu2404i which is 
obvious from the damping of Kph at low T. 
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FIG. 5. Magnon mean free paths of Ca9La5Cu2404i (open 
symbols) and Sri4Cu24 04i (full symbols) as a function of tem- 
perature. 



In contrast to the results at low T we do find a pro- 
nounced doping dependence of the magnon scattering at 
high T. The magnon mean free path is much smaller for 
the doped spin ladders of Sri4Cu2404i and of the order of 
the mean hole distance in the ladders (« 25 A) at 300 K 
|]l8| . In the same crystal resistivity measurements show 
a pronounccd change of their temperature dependence at 
T c± 220 K. This anomaly which is frequently interpreted 
in terms of charge ordering is absent in CagLa5Cu2404i. 
Thus, one might argue that instead of the hole content, 
the hole mobility is crucial for the scattering of magnons, 
in order to explain the weak doping dependence of Imag at 
low T and the strong differences at high T. Yet, holes and 
their mobility are not the only source of temperature de- 
pendent magnon scattering as demonstrated by our data 
for the crystal with undoped spin ladders. In this crystal 
we find a less pronounced, but qualitatively similar de- 
crease of Imag for T > 100 K. At room temperature the 
number of both, thermally excited magnons and phonons 
is already rather large. For example, the inverse density 
of magnons at 300 K as calculated for Aiadder = 430 K 
can be interpreted as the "mean distance" of magnons 
and is already one order of magnitude smaller than Imag- 
Further experimental and theoretical studies are neces- 
sary in order to determine the importance of magnon- 
defect, magnon-phonon and magnon-magnon scattering. 

In summary we have measured the thermal conductiv- 
ity of hole doped and undoped spin ladders realized in the 
compounds Sri4Cu2404i and Ca9La5Cu2404i. In both 
cases we find a huge contribution to k due to magnon heat 
transport. We have applied a simple kinetic approach to 
describe this contributions and have obtained the spin 
gaps of the ladders and the temperature dependence of 
the mean free path from our analysis. Obviously, at low T 
scattering of magnons on holes is rather ineffective since 
in this temperature region Umag of doped ladders is al- 
most identical to Kmag of the undoped ladders. For both 
we find a very large mean free path of magnons of sev- 
eral thousand A. At high T increased mobility of holes 
causes a strong damping of magnon heat transport in the 
hole doped ladders due to increased magnon-hole scatter- 
ing whereas Kmag in the undoped ladders only decreases 
slightly. This is reflected in the mean free path which in 
the former case at 300 K reduces to Imag ~ 60A which 
is in the same order of magnitude than the mean hole 
distance (« 25 A). For undoped ladders Imag ~ 500A at 
300 K. This large value is one order of magnitude larger 
than the inverse density of magnons at this temperature 
and might indicate only weak magnon-magnon scattering 
effects. 

This work was supported in part by the Deutsche 
Forschungsgemeinschaft through SFB 341 and SP 1073 
as well as under Grant No. BR 1084/1-1 and BR 1084/1- 
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